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ABSTRACT 

We present the scientific motivation and observing plan for an upcoming detection survey for debris 
disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) 
survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and 
M spectral classes) to the 850 /im extragalactic confusion limit to search for evidence of submillimeter 
excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 
and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of —40° 
to 80°. In this survey, no star will be rejected based on its inherent properties: binarity, presence 
of planetary companions, spectral type or age. The survey will commence in late 2007 and will be 
executed over 390 hours, reaching 90% completion within two years. This will be the first unbiased 
survey for debris disks since the InfraRed Astronomical Satellite. We expect to detect 125 debris 
disks, including ~ 50 cold disks not detectable in current shorter wavelength surveys. To fully exploit 
the order of magnitude increase in debris disks detected in the submillimeter, a substantial amount of 
complementary data will be required, especially at shorter wavelengths to constrain the temperatures 
and masses of discovered disks. High resolution studies will likely be required to resolve many of the 
disks. Therefore, these systems will be the focus of future observational studies using a variety of 
observatories, including Herschel, ALMA, and JWST, to characterise their physical properties. For 
non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar 
dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for 
the nearest stars in the sample (~ 2 pc). 

Subject headings: stars: circumstellar matter — submillimeter, stars: circumstellar disks 



1. INTRODUCTION 
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Debris disks are the dust disks found around many 
nearby main sequence stars. The dust is short-lived and 
so must be continuously replenished by the destruction 
of comets and asteroids in these systems, deduced to lie 
in fairly narrow belts between 10 to 200 AU from the 
host stars. The InfraRed Astronomical Satellite (IRAS) 
was the first and only large unbiased survey of debris 
disks, showing tha t they o ccur around ~ 15% of n earby 
stars dB ackman & Parescelll993t iPlets fc Vvnckieilll999l : 
lAumann et al.l 11983 ). There is evidence of a substan- 
tial population of disks too cold to have been detected 
by IRAS and w hich are only accessible to submillime- 
ter observations (jLestrade et al.ll2006HWvatt et al.ll2003^ 
surrounding a 5-15% of s tars. This is substantiated by 
the results of iRhee et all ()2007D who note that - 10% of 
stars show 70 /Ltm dust emission in Spitzer MIPS maps, 
but below IRAS sensitivity (either due to low mass or 
cold dust) , suggesting the overall frequency of disks may 
be as high as 25%. Based on an assumed disk incidence 
of 25%, our survey of 500 stars could yield as many as 
125 disk detections, of which « 50 could be cold disks. 
The survey is of sufficient size to measure the disk fre- 



Buildings, The Parade, Cardiff CF24 3AA, U.K. 

Rutherford Laboratory, Space Science & Technology Depart- 
ment, CCLRC Rutherford Appleton Laboratory, Chilton, Didcot, 
Oxfordshire OXll OQX, U.K. 



2 



Matthews et al. 



quency across spectral type, and as a function of age and 
multiplicity. 

Submillimeter observations have already been pivotal 
in studies of debris disks, having imaged or discovered 
seven of the fourteen resolved disks and constrained the 
mass and temperature and hence radial extent of many 
of the remainder. The submillimeter emission is both 
optically-thin and sensitive to relatively large, cold grains 
which dominate the massses of these disks. Hence, our 
survey will provide immediate mass estimates for de- 
tected disks. Constraining the temperature of newly 
identified disks will require complementary data in the 
mid- or far-IR. This will be particularly critical for warm 
disks, since the two potential submillimeter wavelengths 
observable through our survey (850 and 450 fim) will 
both lie longward of the peak of the spectral energy dis- 
tribution (SED) and provide little constraint on the tem- 
perature. 

The study of these disks is revolutionizing our under- 
standing of planet formation. For the fourteen disks 
which have been resolved, observed structures have even 
been used to infer the location of unseen planets (e.g., 
lWvattJ[200l . Many more disks have been character- 
ized by their SEDs, showing that they are the extrasolar 
equivalents of the Kuiper and asteroid belts of the So- 
lar System. The radial positions and masses of these 
belts, particularly when this information can be com- 
pared for stars of different ages, spectral types, mul- 
tiplicity or known planetary companions, provide vital 
constraints on planet formation processes and on how 
the resulting planetary systems su bsequently evolve. 

SCUBA-2 (jHolland et al.ll2006D is a new submiUime- 
ter camera arriving at the James Clerk Maxwell Tele- 
scope (JCMT) in late 2007. SCUBA-2 is expected to 
be, per pixel, five times as sensitive as its predecessor 
SCUBA (S ubmillimeter Common User Bolometer Ar- 
ray, Holla nd et all Il999| ) at 850 /^m. Its larger (7' x 
7'), fully-sampled field of view makes it a premiere sur- 
vey instrument, and seven comprehensive surveys have 
been planned to maximize the scientific output from the 
JCMT over the next five years ( War d-Thompson et al. 
2007, in prep; SLS paper (Plu me et al.ll2007l ). Like its 
predecessor, it observes simultaneously at 850 and 450 
fim, with resolutions of 15 and 7", respectively. Unlike 
its predecessor, SCUBA-2 will Nyquist sample the sky at 
850 fim; at 450 /zm, the array is undersampled by a factor 
of two. The focal plane area coverage is provided by four 
"sub-arrays" at each wavelength, w here each sub- array i s 
comprised of 32 x 40 detectors f see [Holland et al.]|2006D . 

The SCUBA-2 Unbiased Nearby Stars (SUNS) survey 
will utilize 390 hours to observe 500 nearby stars (the 
100 nearest M, K, G, F and A stars) to the JCMT's 
extragalactic confusion limit at 850 /im to detect and 
map circumstellar dust. The survey is completely unbi- 
ased; no star will be rejected due to its intrinsic proper- 
ties. The survey will determine the incidence of disks 
around nearby stars, constrain masses (and tempera- 
tures of disks detected in the far-IR), discover disks too 
cold to be detected in the far-IR, and provide limits on 
the presence of dust which are vital to targeted planet 
search missions such as Darwin and the Terrestrial Planet 
Finder (TPF), as well as future missions which will re- 
solve disks in unprecedented detail, i.e., the Atacama 
Large Milliemeter-submillimeter Array, the James Webb 



Space Telescope, an d the Far-IR Interferometer (FIRI, 
llvison fc Blainll2005l ). 

The mass sensitivity of the survey is a strong func- 
tion of both distance and disk temperature. For disks 
of 70 K, the sensitivity ranges from 0.003 — 1.4Miunar 
(A^iunar = 1/81 Mq) for the nearest (2 pc) and furthest 
(45 pc) stars. For lower disk temperatures of 40 K, the 
range is 0.005 — 2.7Miunar- At the mean distance of the 
survey stars (15 pc), we will be sensitive to dust masses 
typically ~ 200 times the dust mass of the Kuiper Belt 
(« lO^^M®), which is the mass of the disk around e Eri- 
dani. Thus, while it will not be sensitive to present day 
Solar System analogues, it will detect such systems which 
are in a period of unusually high dust mass. We know 
that the Kuiper Belt is a factor of w 100 times less mas- 
sive than expe cted from the dis tribution of solids in the 
Solar System (|Morbidellill2004D . Thus the Kuiper Belt 
used to be more massive. Equally, an episodic mode of 
dust creation could render the Kuiper Belt detectable in 
the future. The mass limit for undetected systems will 
also be useful for future planet search missions. Missions 
such as Darwin/TPF require a dust free system, typi- 
cally below 10 times that of the Solar System, to limit 
the integration time re quired to detect Earth-like planets 
(|Beichman et al.ll2004[) . 

The SUNS survey will provide a significant legacy for 
the field of extrasolar planetary system research. This 
field is rapidly evolving, and a variety of techniques are 
being developed to characterize the planetary systems 
of nearby stars. This survey will determine the dust 
and planetesimal content of these systems and so pro- 
vide vital complementary information on the outcome 
of planet formation in them; for some techniques, the 
dust content even provides the limiting factor determin- 
ing whether such techniques are going to work (e.g., TPF, 
Beichman ct al. 2006b). 

In this paper, we provide details about the SUNS sur- 
vey, including the motivation for the survey in the con- 
text of our current understanding of debris disk sys- 
tems (§ [2|), advantages of the submillimeter compared 
to shorter wavelengths (S 12. 2|) . and the mass sensitivity 
of our survey (§ 12 3p . Our science goals are described 
in § m and we describe the details of our target list in 
§ m including ancillary targets, subset populations and 
statistics. Plans for complementary data to the SUNS 
survey are described in § [5l We describe the data prod- 
ucts in § [6] and summarize the survey in § [T] 

2. MOTIVATION FOR THE SURVEY 

2.1. Current Status of Debris Disk Studies 

All of the approximately 200 known candidate de- 
bris disks were first discovered by their thermal emis- 
sion, which is brighter than the photospheric emission of 
their host stars at far-IR and longer wavelengths (e.g., 
[Mannings &: Barlowlll998h . The majority of these disks 
and disk candidates were discovered by IRAS, which 
provided the first and only large unbiased survey of 
nearby stars for excess thermal emission at 12-100 /Ltm 
(with resolutions 0.5' to 2'). This survey showed that 
~ 15% of nearby stars exhi bit detectable excess emission 
(jBackman fc Parescdll993[ ). and the IRAS disk candi- 
dates have been the subject of intense follow-up obser- 
vations from the ground at a range of wavelengths from 
optical to millimeter. Re-analysis of the IRAS and ISO 
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databases have resuhed in additions, and notably revi- 
sions, t o the stars with id entified infrared excess; for in- 
stance, iRhee et all (I20Q1 found 153 IRAS excess stars 
among Hipparcos dwarfs, 37 of which are newly identi- 
fied. Forty-eight of their excess stars ar e among the 60 
disk candidates identified bv .Moor et alJ (,2006) through 
a re-analysis of IRAS and ISO targets, eleven of which 
were previously unknown. Based on these analyses, Rhee 
et al. revised the fraction of nearby A stars with detected 
60 /im excess to 20%. 

Submillimeter observations have been pivotal in follow- 
up studies of these disks: imaging with SCUBA 
had mapped or discovered seven of the ten resolved 
disks at the time of its d ecommissioning in mid- 
2005 (e.g., lloiiandiEan IMl). Such images have also 
ruled out "disk candidate" stars which turned out to 
have nearby, unassociated background sources which 
fell within the relatively large IRAS beamsizes (e.g., 
iJavawardhana et al.ll2003 ). Finally, photometric submil- 
limeter observations of many disks have provided the best 
contraints on disk masses, as well as constraining the 
SED and theref ore the temperatur e and radial extent of 
the disks, (e.g.. iSheret et al.ll2004[ ). 

Subsequent surveys have been more modest than IRAS 
in the number of stars surveyed but have probed more 
deeply in sensitivity. For example, both the ISO and 
Spitzer strategies targeted several well-defined samples, 
each of which comprised at most 100 stars. The emerg- 
ing picture is that the fraction of star s with detectable 
disks is a f unction of both stellar age (Rhee et al."2007| 
ISpangler e t al. 200 ij), spectral type (feabing et al. 200 H). 
and w avelength (Laureiis et al.l 

moW- IWerner et al.l 
(|2006f ) summarized the early results from Spitzer sur- 
veys and targeted observations, including two Legacy 
Surveys "The Formation and Evolutio n of P la.net ary 
Systems" ( FEPS, Meye r et al. 2004; Kim et all [20051 : 
iStauffer et all 120051: fHines et all I2OO60 and "Galac- 
tic Legacy I nfrared Mid-Plane Survey Extraordinaire" 
(GLIMPSE. lUzpen et aLll2005D . Surveys have also tar- 
geted stellar associations, such as the nearby ~ 10 

Myr old TW H ydrae Association (Uchida et all 120041: 



Low et al.ll2005l l. Cep 0B2 clusters ( Sici ha- Aguilar et al.l 
2000), Orion OBla and OBl b fflernandez et al.ll2006f) . 



and the ope n clusters M47 (iGorlova et all l2004f) an d 
the Pleiades llGorlova et al] 120061 : IStauffer et al.ll2005D . 

Projects to target near by stars have concentrated 



A stars (Rieke et a l. 2005*) and 
(|Brvden et al 



. field F GK stars 

20061 : [Beichman et al.l l 2006aD. includ- 



ing selectio n by known planets (Bcichm an et al.l l2005f ) , 
metal hcity ([Be ichman et al.l [2006b), IRAS 60 um ex- 
cess (iCheneTa l. 2006) and age ([Smith et all I2006D . 
[Beichman et al.l (|2006b( ) concluded that there is no in- 
creased incidence of disks around stars of increased 
metalhcity, in agreement with earlier results based 
on d isks detected by IRAS ([Greaves. Fischer fc Wvatl] 
[2006h . 

Evidence of a disk has been detected ar ound the 
white dwarf at the center of the Helix Nebula ([Su et al.l 
|2007| ). Spitzer has also produced serendipitous de- 
tections of disks around ind ividual stars, including 
HD 4 6190 ([Sloan et al.l I2004f) and Be stars in the 
LMC (Kastne r et al.l I2OO60 . Targeted observations to- 
ward the well-studi ed debris disks around F omalhaut 
([Stauffer et al.ll2005D and Vega (|Su et al.ll2005[ ) have also 



been undertaken. 

In general, few extensive ground-based searches have 
been made for excess emission. Recently, however, sub- 
millimeter surveys unbiased by previous far-IR detec- 
tions have shown that there is a substantial population of 
disks only accessible in the submillimeter , since they are 
too cold to have been de tected by IRAS ([Lestrade et al.l 
[20061: IWvatt et al.l [20031 ). As Table [I] shows, while in- 
dividual submillimeter surveys are typically small (10-20 
stars), the detection rate in each is 5-25%, indicating that 
this population is both real and as numerous as the disks 
which were detected by IRAS. Statistics are currently too 
poor to ascertain whether the presence of such cold disks 
is favored around young/old/early-type/late- type stars. 
[Lestrade et al.[ (|2006) derived a detection rate of 13tg % 
for M dwarfs between 20-200 Myr based on the com bined 
results of their survey and that of [Liu et"aLl ([2004[ ). 

The current debris disk detection rates are shown in 
Figure [1] comparing the detection rates in the far-IR to 
those in the submillimeter. The uncertainties in the disk 
frequencies are still large for several spectral types. Al- 
ready, it is evident that the submillimeter is the favorable 
regime for detection of disks around low-mass stars. 



liir-IR 
(siib)mm 



ipeuiral 
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Fig. 1. — The frequency of debris disks as a function of spec- 
tral type based on existing data. The errorbars are presently 
large for several spectral types. It is also evident that the sub- 
millimeter is the best wavelength range for which to search for 
disks around later-type stars. This is consistent with the detec- 
tion of cold disks. The references for the data are (Greaves et 
aL, in prep, Gauticr ct al.''2007': 'Su ct al.''200g: 'Beichm an et al. ' 
2006a; Brvdcn ct al. 2006; Lestrade ct al. 2006; Beichman et al. 
2005; Naiita & Williamsl 120051 : ILiu'et al.ll2004 IWvatt et al.l 120031 : 
Holmes ct al. 200fK 



2.2. Advantages of the Submillimeter to the Far-IR 

To understand the origin and frequency of debris 
disks, we require an unbiased survey. In the past, 
large surveys for dis ks have been conducted at near- 
IR wavelengths (e.g.. [Mamaiek. Mever fc LiebertI [2002[ : 
[Haisch. Lada fc Ladal[200jl ): however, these only sample 
warm dust in the inner « 1 AU. In the far-IR, maj or sur- 
veys with IRAS, ISO a nd Spitzer (e.g., [S panglcr et al.l 
[20011 : [Decin et al.[ [Werner et alT [20061 and refer- 

ences therein) sample the dust at ~ 50-100 K (typically 
40-100 AU radius), so these are also biased towards warm 
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TABLE 1 

SUBMILLIMETER DEBRIS DiSK SURVEYS 



Reference 


Targets 




-'Vdiaks 


% 


850 fim 1 a rms [mjy] 


Wyatt et al. (2003) 


Lindroos binaries 


22 


3 


13 


1.2 


Holmes et al. (2003) 


Nearby bright stars 


11 


1 


9 


7 


Liu et al. (2004) 


f3 Pic comoving groups 


8 


2 


25 


2 


Carpenter et al. (2004) 


FEPS nearby stars 


127 


4 


3 


28 


Greaves et al. (2005) 


Nearby G Stars 


13 


2 


15 


1.5 


Lestrade et al. (2006) 


M Dwarfs 


32 


2 


6 


0.7 



Note. — Flux limits at 850 fim were adjusted from other wavelengths through scaling of the 
fiducial 850 /jm opacity of 1.7 cm-^ g"'^ and the scaling relations k oc (1/A) and Fi, oc i/(2+/3) 
where /3 = 0.7. 



dust disks. Yet we know cool debris disks do exist; Fig- 
ure [2] shows the SED of the bright nearby disk around e 
Eridani, and Figure [3] shows the fitted temperatures of 
known debris disks versus stellar luminosity for stars of 
differing spectral classes. Disks around lower mass stars 
tend to be significantly cooler than their counterparts 
around A stars, with virtually all such disks exhibiting 
temper atures of < 70 K. This is con sistent w ith the find- 
ings of IWvatt et"al] ([2003il and iRhee et alJ ([2007) of a 
disk population below the sensitivities of IRAS observa- 
tions. 
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Fig. 2. — 850 /im image of e Eridani ll Greaves et al.l2005l) and the 
SEDs of the star (solid line) and grey body fit to the dust (dotted 
line). Note the much greater disk to star contrast in the submil- 
limeter compared to 60 fira, even for this mid-range temperature 
of 55 K. 

To sample the cool dust, we must go to submillimeter 
wavelengths. The number of unbiased surveys conducted 
in the submillimeter, however, is very limited (see Ta- 
ble [T]). In most cases, once the submillimeter detections 
were made, only a few subse quent detailed an alysis of far- 
IR data showed detections (|Zuckermanll2001t ) , indicating 
that the dust is indeed cool. A further point illustrated 
in Tabled] is that the percentage detection rate increases 
sharply with sensitivity, from 3% in the least sensitive 
survey up to 10-15% in deep searches (and higher in the 
P Pic group, perhaps because these stars are young). 
Extrapolating to still deeper surveys, we might expect 



stellar Luminosity [l^(,i„] 

Fig. 3. — The observed temperatures of debris disks versus stellar 
luminosities. The A stars (dark blue) are a sample of 46 detected at 
both 24 and 70 ^lm (W va:tt et al. 2007). The data for the F (light 
blue), G (gre en), K (orange) and M stars (red) are taken from 
IjLestrade et al. 2006, and references therein), with the addition of 
TWA 7 from Matthews et al. (2001). 

a substantially higher detection rate of about 15%. In 
summary, there is growing evidence of a significant pop- 
ulation of submillimeter-bright disks with T < 40 K. The 
final column in Table [T] shows depths reached thus far 
(where equivalent flux limits of other wavelengths have 
been corrected to 850 fim). The SUNS survey, with 
an rms sensitivity of 1 tr = 0.7 mJy will be 40 times 
deeper in flux than the only previous large scale survey 
([Carpenter et al.l[2005f ). This is the equivalent depth to 
the observation of the disk around the 1 Gyr old star r 
Ceti (Figure HI). 

Below, we discuss the effective limits on the submil- 
limeter survey compared to the current generation of far- 
IR detectors on Spitzer. 

2.2.1. Photospheric Contributions 

For the nearby stars where photospheric emission dom- 
inates (Figure [2|) , the dust detection limit of far-IR in- 
struments are set by the accuracy with which a disk ex- 
cess can be discriminated from photospheric contribu- 
tions. The limits on the detectability of a disk are then 
dependent as well on the accuracy of synthetic stellar 
spectra, which are not well-constrained by observations 
at these wavelengths. In this case, it is not possible to 
improve the detection limits, and, in the case of Spitzer, 
detections are limited to > 50% of the photospheric flux 
(|Beichman et al.ll2005l ). As the e Eri SED shows, the 
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Fig. 4. — An example of disk extraction II Greaves et al.ir2004bl) 
is the disk detection toward t Ceti with SCUBA at 850 lira. The 
rms in this map is 0.7 mjy (per smoothed 20" beam). The faintest 
detected peaks are 2 mJy beam~^. The disk is recognized by its 
elliptical shape centered about the star. This star is nearby but 
the disk is in fact one of the smallest yet measured. In addition, 
this field is unusually rich in background sources. 

excess flux over the photosphere is larger by a factor of 
4 at 850 /im compared with 70 /im, and this contrast 
becomes much greater for cooler dust. Note that for the 
very bright AO star Vega at only 8 pc, the photosphere is 
only 5 mJy. Therefore, assuming SCUBA-2 calibration 
to 10%, the largest photospheric error in the survey at 
850 /im will be ±0.5 mJy, so the calibration error will al- 
ways be less than the confusion limit fluctuations of 0.7 
mJy rms. 

2.2.2. Background Confusion Limit 

Background confusion is the absolute flux limit for any 
disk survey. For Spitzer, the confusion limit is 0.7 mJy 
beam~^ at 70 fim (the best wavelength for debris detec- 
tion), but the photometric technique actually limits this 
to « 2 mJy (Beichnian ct al. 2005) and often only 5 mJy 
is reached because of Galactic foregrounds. For 5 mJy 
at 70 /im and grain emissivity of (3 = 0.7, the SUNS sur- 
vey is intrinsically more sensitive to dust cooler than 30 
K. For far-IR limits of 2 mJy, single temperature disks 
would be detectable down to 26 K and for a range of /3 
of 1.0 to 0.0, the limit Hes between 26-43 K. 

2.2.3. Cirrus Confusion 

The 20" beam of Spitzer at 70 /im and the bright and 
complex cirrus background are further limiting factors 
to the far-IR sensitivity. With SCUBA-2, we lose less 
than 20% of stars (i.e., those too far south) and impor- 
tantly, with absolutely no additional bias due the cirrus 
background, which is not a factor in submillimeter ob- 
servations. 

2.2.4. Extended Disks 

Submillimeter single dish observations of debris disks 
are constrained by relatively poor resolution compared 
to optical observations. The beamsize of the SUNSS 
observations will be 15", the JCMT's resolution at 850 
/im. However, the potential does exist to follow-up de- 
tections with 450 fj,m observations with a substantially 
better resolution of 7". Even out to a distance of ap- 
proximately 30 pc, we might expect significant extended 
emission based on the largest examples with 300 AU 



radii (i.e., 20" across) (jSheret et al.ll2004D . and so ac- 
curate photometry requires fully-sampled maps to the 
same limiting sensitivity. Existing surveys are limited by 
the sparse pixel sp acing of current instrumentation (e.g., 
IWvatt et alll2003f ). SCUBA-2's Nyquist sampHng will 
be a valuable new feature, which is important if we are 
looking for relatively cool extended dust around nearby 
stars. 

Our detections of new debris disks will provide some 
clear possible foUowup observations. Filling in the SED 
and thereby refining dust temperatures and spectral in- 
dices can be carried out in the short submillimeter and 
far-IR, and this prospect is imminent through the current 
AKARI all-sky survey and targeted observations with the 
Herschel Space Observatory using its photometric instru- 
ments SP IRE (the Spectral and Photometric Imaging 
REceiver. [G"rifhn et al.l[200l and PACS (the Photode- 
tecto r Array Camera and Spectrometer, iPoglitsch et al.l 
[20061) . 

2.3. Mass Sensitivity 

While the SUNS survey dataset will have a uniform 
flux sensitivity, the mass sensitivity of the survey is a 
function of three factors: the distance of the target, d, 
the temperature (s) of the material in the disks, T^, and 
the opacity of the disks, k^. 

F (f 

At submillimeter wavelengths, in the Rayleigh-Jeans 
limit, the mass of the disk becomes a linear function of 
the temperature in the disk. This makes the mass rel- 
atively straight-forward to estimate if the temperature 
can be constrained. The opacity in the disk is wave- 
length dependent (k^ cx A^^) and difficult to measure 
observationally, but the typical value adopted for debris 
disks at 850 //m is 1.7 cm^ g~^, based on a modified 
black body with /3 = 1 pent et al.l 120001: iPollack et all 

mi. 

Figure [5] shows the mass sensitivity of the SUNS sur- 
vey as a function of temperature for circumstellar dust 
emission at several key distances in our survey. The mass 
of the closest known debris disk, e Eridani, and the mass 
corresponding to ten times that of the Kuiper Belt are 
also indicated. This mass sensitivity is for unresolved 
sources. For resolved sources, our mass sensitivity will 
be lower. 

For new disks, we won't immediately be able to con- 
strain the disk temperature, unless measurements (even 
resulting in non-detections) have been made in the mid- 
IR or far-IR (e.g., through Spitzer, ISO or IRAS). For 
some cold disks (< 50 K), additional 450 /im data will 
provide some temperature constraint, but for warmer 
disks, both submillimeter measurements will lie on the 
Rayleigh-Jeans tail of the SED, and shorter wavelength 
data will be essential. 

3. SCIENCE GOALS 

The outcome of this survey will be an order of mag- 
nitude improvement in the number of disks known in 
the submillimeter. These discoveries will provide a sig- 
nificant legacy for the planet formation community by 
achieving five key goals. The SUNS survey will: 
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Fig. 5. — The disk mass sensitivity (3 cr) as a function of dust 
temperature is shown for various limiting distances. The outer 
bounds of the A stars sample is roughly 40 pc. The limits for F 
and G stars are 25 pc, and all the M stars lie within the 10 pc 
boundary. The horizontal lines show the mass of the e Eridani disk 
(solid) and ten times the mass of the Kuiper Belt (dashed). An 
opacity of 1.7 cm^ g^^ is adopted. 

• determine unbiased statistics on the incidence of 
disks around nearby stars; 

• constrain masses and temperatures of disks previ- 
ously detected in the far-IR (e.g., by IRAS, ISO, 
Spitzer, AKARI, Herschel); 

• discover numerous disks too cold to currently de- 
tect in the mid- and far-IR, for which complemen- 
tary data at shorter wavelengths will be critical to 
constrain temperatures and hence masses; 

• be the basis of source lists for future observing cam- 
paigns using instruments such as ALMA, Herschel 
and JWST; and 

• provide limits on the presence of dust around 
nearby stars that are vital to future missions such 
as TPF. 

Planet formation is thought to be largely complete by 
10 Myr, so studying debris disks tells us where planetes- 
imals remain after planet formation processes are com- 
plete. Such images may indicate locations where planets 
may be located within the systems. The majority of de- 
bris disks we detect will be the extrasolar counterparts to 
the Kuiper Belt. Thus, their temperatures, which will be 
determined through complementary far-IR data, tells us 
the sizes of the planetary systems which are clearing the 
inner regions of dust and planetesimals. Submillimeter 
fluxes also provide the mo st reliable method for measur- 
ing the mass of the disks ()Zuckermanll200ll ). 

We can thus address questions of what affects the ex- 
tent of a planetary system and the mass of its plan- 
etesimal component. For example, do more massive 
stars form the most massive planetary systems? Are 
the disks truncated uniformly with spectral type in- 
dicating that a universal process is at play, such as 
close stellar encounters in the cluster in which the star 
formed, or photo-evaporation of the outer disk edge by 
a nearby massive star in such a cluster. It already ap- 

J iears that more rnassiv e stars host more massive disks 
Greaves fc Wvattll2003f ) that are more readily detected 



(jHabing et al.ll2001h . That said, a comparison of the 
disk masses around members of the nearby young TW 
Hydrae Association suggests no correlation between disk 
mass and spectral typ e for these reasonably coeval disks 
(|Matthews et al.ll2007f) . In addition, comparisons of disk 
mass es around low m ass stars of c omparable age, i.e., AU 
Mic (|Liu et al.ll200l and TWA 7 ([Matthews et al.ll2007f) . 
reveal a difference of a factor of 10. Improved statis- 
tics are critical to reveal whether a genuine trend ex- 
ists. Moreover, information about disk radii is scarce at 
present - the only reliable determinations a re for the few 
disks which have submillimeter detections (jSheret et al.l 
l2004f ). and these comprise a significantly biased dataset. 

With the SUNSS data, it will be possible to test for 
planetary system evolution. Based on our understanding 
of the Solar System's evolution, it is reasonable to expect 
that both the radius and mass of the belts will change 
with time. Within our own Solar System, the outer 
planets may have migrated from the locations where 
they formed (|Malhotralll995l ). a scenario used to explain 
structure in some debris disks ([Wyatt 2003). For exam- 
ple, Uranus and Neptune may have formed in between 
Jupiter and Saturn until their orbits became unstable 
at a few hundred Myr (jThommes. Duncan &: LevisonI 
|1999| ). The radius of the Kuiper Belt may have also 
been significantly different in the past. The period of late 
heavy bombardment in the inner Solar System and other 
evidence, such as the fact that t he Kuiper Belt m ust have 
been more massive in the past ()Morbidellill2004f ) . also in- 
dicates that the dust flux from the Kuiper Belt must have 
varied significantly with time. Thus, understanding the 
evolutionary changes in extrasolar systems will not only 
help place the evolution of the Solar System in context, 
but could also help to illuminate Earth's own origins. 

The survey will test the applicability of the Solar Sys- 
tem itself. One model for the origin of debris disks pro- 
poses that they flare into detectability when a system 
matures to the extent th at a Pluto-sized object ca n form 
in its planetesimal belt (jKenvon fc Bromlevll2004 ). The 
new planet then perturbs the belt out of which it formed, 
initiating the coUisional cascade that produces the co- 
pious amounts of dust that we see for a short amount 
of time. Since such Pluto-analogies may take longer 
to form further from the star, it can take up to 1 Gyr 
for those outer disks to become detectable. This model 
would mean that we could expect to detect disks at larg e 
distances around stars at late times (iRieke et al.ll2005f l. 
Short bursts can also be produced when lunar-sized em- 
bryos are ejected from the inner planetary system. Such 
processes may explain the presence of dust around old 
stars, as well as the app arently episo dic nature of the 
debris disk phenomenon (jRieke et al. ''2005'). Such disks 
would be cool and so may only be detectable at sub- 
millimeter wavelengths. Indeed it can be shown that 
for se veral disks the dust we are seein g must be tran- 
sient ()Su et al.ll2005l : IWvatt et al.ll2007D . An alternative 
view is that disks may be quasi-static in radii. These 
models say that dust mass (but not radius) does evolve 
with t ime, decaying through t he coUisional erosion of the 
belts (jPominik fc Deciiill2003D , and such models can suc- 
cessfully explain the currently available statistics (Wyatt 
et al., submitted). Episodic periods of high dust mass 
are then p rovoked through the destruction of large plan- 
etesimals (|Wvatt fc Dent! [20021 ) or a recent stellar flyby 
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(|Kalas. Deltorn fc Larwoodll2001h . Our survey will con- 
front these models with hard statistics. 

The unbiased nature of the SUNS survey means that it 
will be possible to correlate disk parameters with other 
planet formation indicators. For example, we can as- 
sess whether the presence of a giant planet orbiting near 
the star, or a distant planetary companion detected with 
adaptive optics, affects the mass or size of an outer disk. 
Such planets may form and evolve independently from 
the planetesimals in the system's outer reaches. How- 
ever, evolutionary processes which lead to planets mi- 
grating inward (e.g., massive gas disk, fast core growth) 
may be related to evolutionary processes further out, ei- 
ther promoting or inhibiting the presence of cold disks. 

These are issues which can only be resolved with a large 
unbiased survey so that the variations of disk frequency 
with spectral type, binarity, and age can be untangled. 
The submillimeter is the ideal waveband in which to do 
this, since it is not limited to warm (40-100 K) disks 
(and so small planetary systems). The incidence and 
evolution of mass or radius of cold (< 40 K) disks is 
completely unconstrained at present. Using the submil- 
limeter also means that an estimate for the disk mass is 
obtained immediately, with more accurate limits on tem- 
perature and morphology possible once complementary 
data are obtained at shorter wavelengths. As Figure [6] 
shows, the SUNS survey will be particularly sensitive to 
disks around M stars, since the circumstellar dust is cold 
even at relatively close proximity (< 100 AU) to the host 
stars. 

For the nearest disks, the SCUBA-2 observations will 
result in resolved images. During its lifetime, SCUBA 
imaging at 8 50 /J,m imaged seven of the fourteen re- 
solved disks (iHolland et all Il998t iGreayes et al.l 119981: 
Sheret et al.' '2004"; 'Greaves et alj l2004bl : IWvatt et all 
2005; Williams et al. 2004). Images allow the inner and 
outer radii of the disk to be determined directly, thereby 
confirming inferences about inner cavities. The mor- 
phology of disks c an thereby r eveal the presence of un- 
seen planets (e.g.. IWvattll2003[ ). The orbits, masses and 
even the evolutionary histories of planets have been con- 
strained in this way. Such images will thus be an ex- 
tremely powerful tool for telling us about a region of pa- 
rameter space unreachable with other techniques: plan- 
ets perturbing debris disks are both small (Neptune size) 
and at large distances (tens of AU). 

For disks that are not resolved, the information re- 
turned by SCUBA-2 will provide critical information, 
since it can confirm that the emission is star-centered. 
Associations of excesses with stars are not always re- 
liable from far-IR d ata alone because of the relatively 
large beamsize (e.g.. iJavawardhana et al.r 20021. For un- 
resolved sources, SED fitting and radius estimation will 
indicate those that could be resolved with future obser- 
vatories, at a range of wavelengths from optical to mil- 
limeter. Of particular importance will be ALMA follow- 
up, with our survey giving prior knowledge of the sub- 
millimeter disk flux and good indication of angular size. 
These parameters can only be guessed at from far-IR 
data alone, so with these data far less ALMA time will 
be spent observing sources which are ultimately too faint 
to detect, or performed using non-optimal array config- 
urations. 
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Fig. 6. — These plots show the parameter space in which the 
SCUBA-2 survey is uniquely sensitive {shaded area). In each case, 
this is the location of cold dust. The different lines show the lim- 
its of Spitzer 24 fim sensitivity (dotted line), assuming a limit of 
^disk/-P'* = 0.1 and Spitzer 70 fim sensi tivity (da s hed li ne), assum- 
ing a limit of -Fdisk/^* = 0.55 based on lSu et all H20071 '). The solid 
lines are the SCUBA-2 3 a limits for the closest and furthest stars 
in the samples, assumed to be 2.6-45.2 pc for AO, 3.6-24.8 pc for 
GO, and 1.8-8.6 pc for MO. For AO stars, this space is relatively 
narrow (shown by the shaded region). For M stars, however, the re- 
gion is much more extensive and includes significantly smaller disks 
below the detection limits of Spitzer at 70 fj,m. This figure does 
not account for the fact that at these distances the largest disks 
are going to be significantly resolved, and therefore, we won't reach 
the 2 mjy limit, i.e., the SCUBA-2 limits are optimistic for disks 
larger than the beam size. 

IWvatt et al.l (|2007f ) modeled the population of disks 
aroun d A stars, fitti ng to the Spi tzer statistics of 
Rieke eiaU (|2005[ ) and lSu et all (|2007[ ). yielding predic- 
tions for the outcome of the A star component of SUNSS, 
excluding the cold disks which would not have been de- 
tected by Spitzer. They predict that 17 (of 100 stars) 
would have detected disks at > 2 mJy and that > 5 
would be resolvable in 450 fim imaging. This study also 
concluded that the SUNSS would be vital for constrain- 
ing the distribution of disk radii and for determining the 
unknown population of disks that are too cold to detect 
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with Spitzer (not included in the 17% of expected disks). 
Assuming the other spectral types have a similar sub- 
millimeter detectability, we could justify a prediction for 
the survey of 85 detections and 25 resolved disks based 
on the model of the A star population. A cold disk pop- 
ulation not currently sampled by the Spitzer statistics 
would result in an elevated number of disk detections. 
We have estimated this value as approximately 10% (50 
disks) , but deducing the true incidence of cold disks is in 
fact one of the principle drivers for this survey. 

4. TARGET DETAILS 

4.1. The Sample 

The survey will cover 500 nearest stellar systems - 100 
primaries of each of the spectral types A, F, G, K, and 
M observable from the JCMT. Subgiants are included 
for spectral types A, F and G, while only dwarfs are in- 
cluded in the K and M targets. The aim is to obtain 
samples that are statistically robust and can be inter- 
compared, while keeping the survey completely unbiased 
with regard to choice of star. This is a unique feature of 
our survey; no star will be rejected because of its intrin- 
sic properties. The nature of the mass function of stars 
means that the five sub-samples cover different volumes. 
The distance limits extend out to 45, 25, 22, 16.5 and 
8.6 pc for A, F, G, K and M stars respectively. The un- 
certainties on the parallax for these stars are ±1 pc for 
A, F and G stars and ±0.1 pc for K and M stars. These 
volume explored by this survey are similar to those be- 
ing explored with Spitzer and earlier with ISO, but the 
SCUBA-2 survey will be the first unbiased study. Fur- 
ther details on the selection of targets will be presented 
in a forthcoming paper (Phillips et al., in preparation). 

The allocation of 390 hours includes 330 hours dur- 
ing the first two years of SCUBA-2 operations. Of this 
time, 270 hours has been allocated to weather band 3 
(0.08 < T225GHZ < 0.12) and 60 hours is available within 
band 2 (0.05 < r225GHz < 0.08) weather. (The lowest- 
elevation stars are prioritised for the better conditions.) 
The number of targets that need close observing con- 
straints (having low elevation and thus priority for up- 
per band 2 conditions) is only 10% of the total number of 
targets. The remaining 60 hours is allocated beyond the 
first two years of operation in band 2 conditions; some 
foUowup may be possible in this period, although com- 
pletion of the 850 /im survey will have priority. 

Figure [7] shows the targets are distributed across the 
sky; the furthest distance is 45 ± 1 pc, and there is no 
clustering of stars towards the Galactic Plane. Observa- 
tions at any time of the year will have suitable targets. 
The declination limits of —40° to -1-80° ensure that tar- 
gets rise to at least 30° elevation. 

We will image the total sample of 500 stars down to our 
adopted JCMT confusion limit of 0.7 mJy at 850 /im. In 
borderline band 2/3 conditions (t225ghz = 0.08) at air- 
mass 1.15 (30° from zenith), this sensitivity requires an 
average of 45 minutes observing per star. The priority 
will be for 850 ^m observations with 450 fXTn imaging a 
bonus (giving first hints of disk sub-structure suitable for 
follow-up). In general, 450 ^m observations will form a 
follow-up Pl-driven project once the population of suit- 
able detected disks is known. The practical strategy is to 
start with the nearest targets and work outwards, while 
maintaining balance among the spectral types. Thus the 
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Fig. 7. — Distribution of target sources in the DDS. The targets 
are all-sky. The declination limits are —40 < 5 < -1-80. 

best resolved disks will be observed first, and these will 
have the highest priority for follow-up imaging. 

Figure [5] shows the variation of mass sensitivity for 
disks of different temperatures and the limitations for 
several key distances in the survey, including the inner 
bound (2 pc) and the outer bound (45 pc) , as well as the 
outer bounds of the M dwarf (8.6 pc) and G and F stars 
(22 pc and 25 pc). This figure shows that, for unresolved 
disks, our sensitivity is sufficient to detect masses com- 
parable to that of e Eridani (0.002 Mg, = 0.16 M^nar) 
at most temperatures within 10 pc. The mass of e Eri- 
dani lies in the mid-range of d ust masses for the h andful 
of detected Sun-hke systems ([Greaves et al.|[200^ . The 
detectable mass is of course modified by dust location, 
stellar heating, and resolution, but we note that 15 pc 
is a useful reference point as it is the far end of planet 
search distances for TPF. At that distance, we will de- 
tect all unresolved systems with temperatures exceeding 
70 K. Such limits mean that non-detections result in con- 
straints on dust mass in these systems of typically 200 
times the dust content of the solar system, but down to 
only approximately four times this for our nearest target, 
GJ 699, at 1.8 pc. 

4.2. Ancillary Targets 

In addition to the 500 primary targets of the survey, 
the large FOV of SCUBA-2 will enable us to image com- 
panions in multiple systems, and in some fields, to in- 
clude other nearby stars where they are present. In our 
target list, at least 240 stars are members of multiple 
systems of at least two members. This high fraction of 
multiple systems raises the total number of stars sur- 
veyed to well over 750. Very few of these systems require 
more than one SCUBA-2 sub-array, since even wide bi- 
naries up to 1000 AU in separation will have an angular 
separation of < 2'. In the case of binaries, their inclusion 
in the survey is based on the spectral type of the primary; 
exempting the secondary from the initial sample ensures 
that no bias is introduced due to the probabilities of disks 
in binary systems or the weighting of secondaries toward 
lower masses. 
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4.3. Subset Populations 

The sample-size for each of the spectral types was cho- 
sen so that we could statistically distinguish between de- 
tection rates of 5, 10, 25 and 50% when the primary 
dataset is divided into the following sub-categories. 

• Stellar type, with 100 stars each of A, F, 
G, K and M. (Properties here and subse- 
quently are from t he NASA NStars database at 
http://nstars.arc.nasa.gov[ ). 

• Stellar age, with ^ 150 stars younger than 1 Gyr 
and ^ 350 stars aged 1-10 Gyr. The division cor- 
responds approximately to the end of the heavy 
bombardment phase of the Earth, and this age split 
can readily be applied from the decline in X-ray ac- 
tivity, as measured by the ROSAT all-sky survey 
(jGaidos et all [1991 . The 150:350 division arises 
from the proportions within each stellar type falling 
with this age bracket. This is ^100% of A-stars (1 
Gyr on the main-sequence) , 20% of F-stars (lasting 
5 Gyr) and 10% of G, K and M stars (hmited by 
the Galactic Disk age of 10 Gyr). 

• Stellar multiplicity, with w 50% of target stars hav- 
ing another star as a companion. 

• Presence of a planetary system, with ^ 20 exam- 
ples within our volume known at the present time 
from radial velocity searches. 

4.4. Statistics 

The overall choice of sample size is driven by the aim 
of differentiating detection rates among the listed sub- 
groups. The key is to have a sample size large enough 
that the smallest sub-group with the lowest detection 
rate can be distinguished from rates in other sub-groups. 
The detection rate for X disks among Y stars is defined 
sls Z — {X + / — ^/X)/Y as the uncertainty is Pois- 
sonian for small counts. In the simplest interpretation, 
we accept rates as differing if the upper bound of one 
Z value is below the lower bound of the comparison Z 
value. In practice this means the survey size is driven by 
the smallest samples, that are one-tenth of the compar- 
ison group. The two cases are the planetary systems (~ 
1 in 10 Solar-type systems) and young Sun-like stars (20 
among 200 of the G and K stars). Here we can distin- 
guish observed detection rates Z for any pair of intrinsic 
detection rates among 10%, 25% and 50%. The formal 
minimum required to do so is 163 stars, which is just 
met by our smallest group: the approximately 180 stars 
of late-F to mid-K potential planet-hosts. 

Detection rates do not need to be compared within 
every spectral type. For example, the effects of age on 
debris will be investigated for stars of Sun- like mass, i.e., 
types both G and K, rather than G and K separately. 
In fact, this approach is required to allow us to establish 
robust distinctions (i.e., at the 2-3 a level) between de- 
tection rates as a function of spectral type and age, due 
to the practical limits on sample size. For instance, with 
200 stars per sample, a 15% difference in detection rate 
can be identified at the 3 cr level. 

Results from the far-IR surveys with IRAS, ISO and 
Spitzer show that these rates actually occur amongst 



sub-groups of stars. For example, 10% is th e global 
rate for Sun-hke stars (jGreaves fc Wvattll2003f ). 25% is 
the rate among Sun-li ke stars with known giant planets 
(|Beichman et al.ll2005l). and 50% is the rate among more 
luminous A stars 1g reaves fc Wvatd [200l iRieke et all 
l2005f l. It is useful to note that 10% of disks exist even 
at the limited IRAS sensitivity, and so this a minimum 
frequency. However, we can also distinguish a 5% disk 
occurrence should it occur among one of the less well- 
studied sub-groups. In the smallest sub-groups of only 
20 stars, we have evidence already for detection rates o f 
20-25% ([G reaves fc Wvattll200l iBeichman et al.ll2005[ ). 

5. COMPLEMENTARY DATA 

To make the most of the JCMT Legacy survey dataset, 
it is essential that we ensure that the targets are selected 
based on accurate spectral classification. Although the 
targets are all very near the Sun, this proximity does 
not guarantee that the optical spectrum of each is well- 
characterized. This is critical to exclude K-giants and 
brown dwarfs, as well as to ensure that the sub-categories 
of the targets are accurately populated. A s part of the 
ongoing RECONS survey of nearby stars, iHenrv et al.l 
(2006) have recently published spectral classifications for 
most of the nearby stars on our target list. For those 
which have not yet been observed spectroscopically, we 
have undertaken a small optical survey. We are using the 
1.2 m telescope at the Dominion Astrophysical Observa- 
tory to obtain the bulk of these spectra with a resolution 
R = SX/X w 3000. 

The SUNS survey is complementary to Spitzer pro- 
grams searching for warmer debris disks at 24, 70 and 
160 microns. Synergy with Spitzer data is important 
to construct the spectral energy distributions of the full 
disk population, even where Spitzer has failed to detect 
a disk in our targets. Due to the timescales of the mis- 
sions, most far-IR complementary data for our survey 
targets are likely to come from AKARI and the Herschel 
Space Observatory. AKARI is currently doing an all-sky 
far-IR survey. A targeted Herschel debris disk survey is 
planned using guaranteed time on the Herschel photo- 
metric instruments PACS and SPIRE. In addition, other 
photometric observations are likely to be proposed dur- 
ing open time. 

6. DATA PRODUCTS 

The data products of the survey will include an archive 
of 500 single fields centered on the 500 target stars, pub- 
lication quality images of all the detected disks, and a 
catalogue of the measured fluxes and flux limits at 850 
and 450 /im. In addition, plots of the fitted SEDs for 
each disk and star combination will be made using opti- 
cal and IR data combined with the submillimeter fiuxes. 

It will be possibl e to extract systein atically parameters 
from the data fsee lSheret et al.ll2004f ). These will also be 
included in the Legacy database of the survey, and in- 
clude dust temperatures, masses, spectral indices, and 
characteristic disk radii (either inferred from tempera- 
ture and spectral index, or measured where the disks are 
resolved). Systematic modeling will also yield masses of 
colliding bodies within the debris belts (up to specified 
sizes), and estimates of mass, location and orbital direc- 
tion for perturbing planets (where the disk structure is 
well- resolved) . 
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We have already developed additional tools for mea- 
suring disk sizes and fitting spectral energy distributions 
and will make these publicly available. For examples of 
development of public databases that we have already 
made available^^. Further, our 5piteer contributors have 
developed tools using HIRES techniques, specifically for 
resolution-enhanced mapping of debris disks. These tools 
will be used in a unified analysis of the SGlJBA-2 / Spitzer 
database. Herschel data will be incorporated to our anal- 
ysis where possible. 

6.1. Ancillary Data Products 

Some important ancillary data products will emerge 
from the survey. These will include a catalogue of fluxes 
and positions of background sources not associated with 
the star, the majority of which will be distant, dusty 
;alaxi es of the kind discovered by ISmail. Ivison fc BlainI 
1997| ). and a list of fluxes for nearby stars where the 
photosphere is detected but there is no debris disk. 

The extragalactic catalogue has value in providing a 
large database of submillimeter galaxies (SMGs) that 
can be followed up with adaptive optics to measure op- 
tical/IR properties at very high spatial resolution. A 
unique aspect of these target fields is that the survey 
stars themselves will be available for guiding, and they 
are technically ideal, with magnitudes down to about 
my = 12 and offsets of up to about 5 arcminutes. Fur- 
ther, the fields add up to 7 square degrees of sky ob- 
served to the JCMT's 850 confusion limit but split 
500 ways rather than in a few fields. Thus the sources 
counts provide a test for cosmological variance, in com- 
parison to the SCUBA-2 cosmology survey. Each 49 
square arcminute survey fleld is well-matched to the fleld 
of view of the KMOS near-IR facility on ESO's 8-m VLT 
(jSharples et al.ll2004f ). allowing up to 24 SMGs to be 
targeted spectroscopically using deployable integral fleld 
units, which promises to yield the least-biased spectro- 
scopic redshift distribution to date. 

The stellar photospheric data are valuable because the 
long-wavelength end of the Rayleigh- Jeans tail has never 
been measured. As an example, the AO star Vega at 8 
pc has a photosphere of about 5 mJy at 850 fj,m, and 
so the equivalent for a non-debris star can be detected 
with our survey at the 7a level. The science value is 
in observing the contribution of the photosphere to the 
SED of the star, a quantity that is poorly measured in 
the far-IR, especially for cooler stars (K and M) with 
deep molecular absorption bands. Measuring signals be- 
low the expected levels would show that the photosphere 
does not simply extrapolate as a blackbody, constraining 
stellar atmosphere models. At present, about 5% of stars 
are known to have far-IR sig nals 3 a or more be low the 
photospheric prediction (e.g.. iHabing et al]l2001h . 

7. SUMMARY 

Many stars are surrounded by dusty cold debris disks. 
These are fed by asteroids and comets orbiting the stars. 
Studying the location, mass and morphology of these 
disks provides crucial information about the outcome of 
planet formation in these systems. For the fourteen disks 
which have been resolved at present, observed structures 
have even been used to suggest the location of unseen 



planets. Hundreds more stars have had their disks char- 
acterised from their SEDs showing that these disks are 
the extrasolar equivalents of the Kuiper and asteroid 
belts in the Solar System. 

In this Legacy survey, we will use 390 hours of SCUBA- 
2 time on the JCMT to observe 500 nearby main- 
sequence stars to search for debris disk signatures. This 
survey will be the flrst unbiased one since IRAS, as pre- 
vious far-IR surveys have had to omit many stars. The 
crucial value of the submillimeter is that the stellar pho- 
tospheric signal is irrelevant and so any star can be ex- 
amined. The output of our survey will be robust statis- 
tics on the incidence of debris disks plus discovery of the 
underlying causes (in terms of the stellar environment 
and history). The nearer systems may also be resolved, 
contributing to planetary detection and planning for mis- 
sions such as TPF/Darwin. 

The data products will be unique, comprising deep and 
uniform searches for debris without any bias towards par- 
ticular stellar properties. This has never been done at 
any wavelength, and particularly not in the submillime- 
ter where a new cold population of disks is barely ex- 
plored. The SUNSS wi ll exceed the modest, un biased G- 
dwarf SCUBA survey (|Greaves fc Wvattl[200l by forty- 
fold in stellar numbers while being substantially deeper. 
The SCUBA-2 sensitivity will approach the Kuiper belt 
dust level for the closest Solar analogues; a disk around 
these targets could actually be detected in our survey 
before the equivalent has been mapped around the Sun. 
The survey can never be done better until large far-IR 
telescopes fly in space ~ resolving the disk spatially from 
the stellar photosphere - a prospect considerably down- 
stream of JWST. 

The Science Legacy lies in answers to the flve key out- 
comes: 

• determining unbiased statistics on the incidence of 
disks around nearby stars; 

• constraining masses and temperatures of disks pre- 
viously detected in the mid- and far-IR (e.g., in 
Spitzer surveys) ; 

• discovering numerous disks too cold to detect in 
the far-IR; 

• being the basis of source lists for future observing 
campaigns using, e.g., ALMA and JWST; and 

• providing limits on the presence of dust that are 
vital to future missions such as Darwin/TPF. 

With these answers in hand, we will be able to un- 
derstand for the flrst time the relation of debris disks - 
tracing planetesimals up to tens of km across in orbits at 
tens of AU - to the inner planetary systems detected by 
other methods. The results, especially when combined 
with shorter wavelength data to constrain temperature 
and mass, will test models of planet formation spanning 
across the scale of our Solar System (from inside Mer- 
cury's orbit to beyond Neptune's). The images of disks 
will be followed in the next decade by high-resolution 
imaging that may indicate perturbing planets, even fol- 
lowing their orbital perturbations in real time. The re- 



sults will be vital for the detection of extrasolar Earths 

1® see |http://www.roe. ac.uk/ukatc/research/topics/dust/indentificatiowit lllntOronagraphs. 
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